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Abstract 

High-resolution numerical experiments, described in this work, show that velocity 
fluctuations governed by the one-dimensional Burgers equation driven by a white-in- 
time random noise with the spectrum \f{k)f oc k~^ exhibit a biscaling behavior: All 
moments of velocity differences SnKsO = \u(x + r) — u(x)J = |Attp oc while 

Sn> 3 {r) oc with Ri 1 for real n > 0 (Chekhlov and Yakhot, Phys. Rev. E 51, 
R2739, 1995). The probability density function, which is dominated by coherent shocks 
in the interval An < 0, is V(Au,r) oc (An)“'^ with g ~ 4. A phenomenological theory 
describing the experimental findings is presented. 
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PACS number(s): 47.27. Gs. 



Our recent study 



of the one-dimensional Burgers equation 
du 

'm 



du „ d'^u 
u— = f + u- 



dx dx^ 

driven by a white-in-time random force dehned by the correlation function 



( 1 ) 



f{k, t)f{k\ f) oc D{k) 6{k + k') 6{t - t') (2) 

with D{k) = 6o and Cq = 0(1) was motivated by an interest in dynamical processes 
which involve an interplay between chaotic and coherent phenomena. It has been shown 
that the velocity held u{x,t) consists of random-in-time and random-in-space huctuations 
superimposed on the relatively strong and long-living shocks. Numerical simulations yielded 
the energy spectrum E{k) oc \u{k)f oc k~^ with x = 5/3±0.02, characteristic of Kolmogorov 
turbulence and the Eulerian correlation function C{k,u) = \u{k,u})f oc <h(a;//c^) 
with the dynamic exponent ^ = 2/3. This result shows that in this system the kinematic 
transport of the small-scale velocity huctuations by the large-scale structures is very weak. 
Investigation of the velocity structure functions S'„(r) = [u{x + r) — u{x)]^ = (Aw)"' with 
integer n revealed strong deviations from the Kolmogorov picture of turbulence: all moments 
Sn> 3 {r) oc rG with 1, characteristic of strong shocks. Thus, the system governed 

by (1) — (2) shows both “normal” (Kolmogorov) and anomalous scalings with the latter 
dominated by the coherent structures (shocks). In this work we are interested in the details 
of the probability density functions (PDF’s) characterizing the huctuations generated by 
(1) — (2) and in the role the structures play in the determination of the PDF’s shape. 
The PDF V{Au,r) is dehned such that V{X,r) dX is the probability of hnding a velocity 
diherence Au = u{x + r) —u{x) within the interval (X, X + dX) for inhnitesimally small dX. 
A spectral code with 12288 Fourier modes was used in the numerical experiment. Equation 
(1) with a hyperviscous (instead of viscous) dissipation term was solved. The details of the 
numerical procedure are reported in 

The most prominent feature of Burgers equation is a tendency to create shocks and, 
consequently, to increase the negative velocity diherences Au < 0 and to decrease the positive 
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ones Am >0 0. Thus, strong asymmetry of the curve V{Au,r) is expected. The two- 
point PDF V{Au,r) was measured for a set of separations r covering a variety of scales in 
the system in the following way. The range of variation of the velocity difference, —5 < 
Au/urms < 5, was divided into 10^ bins. The data were collected during a time longer than 
10 largest eddy turnover times (corresponding to O(IO^) time-steps) and were distributed 
among the appropriate bins to generate a histogram. Fig. 1 presents V{Au, r) for the inertial 
range separations r/dx = 200, 250, 300, 350, 400, where dx = L/12288 is the mesh size and 
L = 27 t is the system size. It follows from (1) — (2) that: {Au)^ oc eorlog rkd and that is 
why this PDF has a shifted maximum, approximately at 0 = {Au)/R^^^ ~ 0.5. Here the 
function R{r) dehned as R{r) = J[f{x + y) — f{x)]^ dy, was also directly measured. It is 
fully force-dependent and in a system with viscosity it may be analytically calculated for the 
inertial-range values of separation r, giving R{r) oc (Am)^ oc eor log(rkd), where kd — >■ -|-oo 
is a dissipative cutoff wavenumber. Technically speaking, the system considered in this work 
does not have a real “inertial range” since the mean dissipation rate e = 0{eolog{L kd)) 
depends on the ultra-violet cut-off kd. This dependence, however, is weak and in what 
follows we take kd = 0{{eo/v‘^)^^^). 

The tail of the PDF V{Au,r) for Au < 0 is shown on Fig. 2 for various magnitudes of 
the displacement r in the universal range. One may observe from this hgure that the PDF 
for Au < 0 may be well approximated as: 

'P(Am, r) oc (Am)”'^, (3) 

with q ca 4. We have also found that V{Au,r) for Au > 0 is well htted by the exponential: 

'P(Am, r) oc (4) 

with the constant a to be determined from the theory. The dynamic argument leading to (3) 
will be presented below. The results shown on Figs. 1, 2 are highly nontrivial because the 
observed algebraic decay of the PDF V{Au, r) as Au/ {Au)rms ~oo leads to the divergence 

of the moments S'„(r) for n > 3 for the inviscid case. However, as we also observed, the 
single point PDF V{u) is a very rapidly decreasing function which is close to the Gaussian 
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and that is why the occurrence of shocks with an amplitude Au > [/q ~ {e Llog{L kd)Y^^ is 



highly improbable and one can expect the PDF V{Au, r) to decrease sharply for Au < —Uq. 
This is sufficient for the existence of all moments Sn{r). A full analytical theory leading 
to an expression for V{Au,r), which unifies both asymptotics (3) and (4) will be published 
elsewhere |^. 

To develop a phenomenological theory we assume that the flow can be represented as a 
superposition of coherent and random components. The coherent contribution is visualized 
as a “gas of shocks” and a single structure (shock) can be approximated by the exact tanh- 
solution of the unforced problem 0 . In particular, let us assume that solution for the normal 
(not the hyper-) viscosity case has the form 



The first contribution to the right side of (5) describes the slowly varying coherent “gas of 
shocks”, whereas the second represents the effects unaccounted for by the first term. Here 



respectively. The physical picture behind this representation is the following: the forcing 
produces the low energy excitations which coagulate into ever stronger well separated shocks 



the shock width /j n/t/j, which means that the stronger the shock, the more narrow it is. 

Statistics of the dissipation rate fluctuations were investigated in detail in Ref. 1. It has 
been shown that the energy dissipation takes place mainly 99%) inside the well separated 
strong shocks. Thus, it follows from (5) that the dissipation rate in interval of length r is 



where we denote Yi = {x — a*) Ui/{2u). The principle contribution to the sum comes from 
the strong and narrow shocks, and, therefore, we can neglect the nondiagonal terms with 
i % j. Taking the integral for inertial-range values of r we have 




( 5 ) 



Oj and Ui denote the coordinates of the centers of the shocks and the shock amplitudes 



due to the non-linear interactions. It will be clear below that the detailed shape of the shock 
assumed in (5) is unimportant. The most essential feature of the ta nh -, solution (5) is that 





( 7 ) 
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( 8 ) 



On the other hand, it can be directly shown from (1) — (2) that 

'LUo 



, ( LUo\ 
6r = eoln \-^j ■ 



Introducing the PDF V{U,r) to hnd a shock with amplitude U in the interval of the length 
r we obtain from the last two relations 

rUo . . . /LUo^ 



/ U^V{U,r)dU oceorln 
from which we readily establish the form of V{U,r) 



V 



( 9 ) 



( 10 ) 

Since V{U,r) = V{U) r/L, the relation (10) establishes the shape of PDF V{U) oc to 
hnd a shock of the amplitude U. Note that r/L is the probability to hnd a shock center 
within the interval of the length r which is in turn placed in the larger interval of the length L. 
The low integration limit U ^ u/L, corresponds to the amplitude of the “weakest structure”, 
contributing to e^. Formula (10) is a consequence of relations (7) and (8), and is valid in the 
logarithmic case when the forcing function is dehned by (2). It is only in this case that we 
can establish the form of the PDF. 

Thus, according to the data presented in Fig. 1 and the theoretical considerations devel- 
oped above, the PDF of velocity differences can be represented as: 

V{Au,r) =aR--^ F{^) (11) 

Ra 

in the interval 0(— 1) < x = Au/R^^^ < oo, where function R(r) is dehned above, a is a 
numerical constant and F{x) is a scaling function (see, P) which is assumed to go rapidly 
to zero when |x| is large. In the interval x -C —1, and |Am| < 0{Uq), where the PDF is 
dominated by the well-separated shocks, we have: 



P(A«,r)=6-^, (12) 

(Am)^ 

where 6 is a constant. When Au -C —Uq, the PDF is a rapidly decreasing function of Au/Uq. 
The moments of velocity difference are evaluated readily with the result: 

Sn{r) = [ {Auf V{Au,r)dAu = hnreo— hS„(eoi?)^, (13) 

J Tl o 
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where the amplitudes Bn depend on the shape of the scaling function F{x). The constants 
bn oc (—1)"' for integer n and for the noninteger values of n the structure functions Sn = | Au|" 
so that relation (13) should include the absolute value of the hrst term in the right side. It 
follows from (11) — (13) that all moments S'^(r) with n > 3 are completely determined by 
the upper cut-off in (13) 



Tjn-3 

Sn{r) oc eor — ^ 
n — 3 



(14) 



which is in excellent quantitative agreement with |Q, whereas for 0 < n < 3: 



n 

Sn{r) oc {eoR)^ , 



(15) 



as in the Kolmogorov theory of turbulence @]. Thus, the anomalous scaling of the velocity 
structure functions S'„(r) appears only for n > 3. It should be stressed that, in accord with 
(13), in the logarithmic case considered in this work the contribution from the shocks to the 
moments Sno is smaller than the one from the scaling component of the PDF only by factor 
1 / log(r kd) which makes the experimental investigation of the details of the crossover very 
difficult. 

The prediction (13) has been tested in [|I[. It has been shown that S' 2 n(r) oc ^ 2 n ~ 

0.91 for n > 2, indicating that these correlation functions are dominated by coherent shocks. 
The results of the measurements of the structure functions S'n(r) with n = 1/3, 2/3, . . . , 6/3, 
presented on Fig. 3, are in good agreement with the scaling law (15). The general structure 
of the moments of velocity differences given by expression (13) — (15) is similar to the outcome 
of the recent theories of the random- force- driven Burgers equation by Polyakov [^] (Id) 
and Bouchaud [§] {d — > cxo) cases. The Polyakov theory conhrmed our qualitative argument 
leading to the algebraic decrease of the PDF in the interval Au < 0. 

Fig. 4 presents the PDF of the shock amplitudes. The problem of the shock location was 
solved in the following simple but reliable way. At each spatial point x the local gradient 
of the solution u{x) was measured. Then, if u'{x) > 0, it was assumed that this point x is 
outside of a shock, otherwise x lies inside of a shock. Once inside a shock, one can march in x 
until the gradient becomes zero, and thus the boundaries of the shock may be located, and so 
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forth. Note that the shock amplitude obtained in this way has been corrected to exclude the 
Gibbs phenomenon typical in spectral approximations of discontinuous functions. To reduce 
the statistical noise in V{U) in Fig. 4, a simple smoothing procedure was applied: V{U) was 
averaged over eight surrounding points. The result presented in Fig. 4 demonstrates that 

V{U) oc U~^ (16) 

is observed for all \U /Urms\ > 0.5. The fact that V{Au) ^ V{U) when Am < 0 tells us 
that in this range V{Au) is dominated by the well-separated shocks. This conhrms the main 
assumption of the phenomenological theory presented above. It follows from Figs. 1 — 4 
and relation (12) that the anomaly in the high-order moments results only from the slow 
(algebraic) decrease of the PDF in the interval Au < 0. As was pointed out above, in this 
case one expects a cut-off at some Au ~ Uq- 

We have also investigated the problem (1), (2) driven by the white-in-time random forces 
with D{k) ^ 0 only for k < 5 and D{k) oc P]. The outcome of the simulations in 

both cases revealed the algebraically decreasing V{Au, r) oc r/ IAmI*^ for Au/{Au)rms ^ ~ 1, 
with the exponent q, related to the functional form of D[k). The former case of the large- 
scale driven Burgers equation was investigated in a recent paper by Bouchaud et. al. P 
using a replica trie in the limit of the space dimensionality d ^ oo. Although the scaling 
of the moments of velocity differences, obtained in Ref. 5 is the same as the one observed 
in our simulations, the shape of the PDF in the Id-case, numerically found by us, differs 
dramatically from V{Au,r) = (1 — r) 6{Au — r) + (3 F{Au/Uq), derived in Ref. 5. Here 
F{x) is a scaling function and /? is a number. This means that the physical mechanisms, 
responsible for the anomalous scaling in the one and multi-dimensional systems are different 
and understanding of the transition between the two behaviors is an extremely interesting 
challenge. The detailed theoretical and numerical investigations of the different cases of 
forcing functions will be published elsewhere |P]. 
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Figure 2: The tail of the two-point PDF V{Au,r) (points) for separations r/dx = 

150,200,250,300,350,400,450 within the universal range, plotted on a logarithmic- 
logarithmic scale. The slope of the solid lines is equal to —4. The graphs for different 
values of r are arbitrarily shifted along the vertical axis for clarity. 
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Figure 3: Velocity structure functions | Au|"' for noninteger values n = 1/3, 2/3, . . . , 6/3 (dot- 
ted curves) . Slopes of the linear least square fits (solid lines) from top to bottom respectively 
are: 0.111, 0.222, 0.330, 0.433, 0.531, 0.620. 





